function [ bbox, area ] = boundingBox( points )
    minx = 2^10;
    maxx = 1;
    miny = 2^10;
    maxy = 1;
    for k=1:length(points)
        x = points(k, 2);
        if (x < minx); minx = x; end;
        if (x > maxx); maxx = x; end;
        y = points(k, 1);
        if (y < miny); miny = y; end;
        if (y > maxy); maxy = y; end;
    end
    maxx = maxx + 1;
    maxy = maxy + 1;
    bbox = [minx, miny, maxx, maxy];
    area = (maxx - minx) * (maxy - miny);
